Delay Curves for Calls Served at Random 

By JOHN RIORDAN 

(Manuscript received March 11, 1952) 

This paper presents curves and tables for the probability of delay of calls 
served by a simple trunk group with assignment of delayed calls to the 
trunks at random and with -pure chance call input. These are contrasted 
with the classic results of Erlang {"Erlang C") which are based on service 
in order of arrival. Trunk holding times for both have an exponential dis- 
tribution. The theoretical development for computation of the curves is di- 
rected to the determination of the moments, which seem to be a naiural means 
of simplification. 

1. INTRODUCTION 

One of the classic results in the study of telephone traffic is the for- 
mula for delay given by the Danish engineer A. K. Erlang^ in 1917. This 
is for random call input to a fully accessible simple trunk group with 
the trunk holding time exponential and calls served in the order of arrival. 
A proof for this formula and a set of curves for its use have been given 
by E. C. Molina.2 

In many switching systems it is not feasible to fully realize this ethical 
ideal of first come, first served, and it has long been of interest to de- 
termine delays on another basis. The contrasting assumption is of calls 
picked at raiidom, which is again an idealization but in large offices ap- 
pears to be called for, as a bound for the service actually given. 

The first attempt to formulate the last seems to be that of J. W. 
Meilor.^ While his basic formulation is incomplete, it offers a useful 
approximation to the complete results, particularly in the most interest- 
ing region of heavy traffic, and will be referred to here as the "Mellor 
approximation." A complete formulation due to E. Vaulot"* appeared 
in 1946 and included both the fundamental differential recurrence rela- 
tion and formulas for delay probabihties for small delays. For complete- 
ness, these are repeated below. F. Pollaczek^ has given a development of 
Vaulot's work directed toward determining an asymptotic delay for- 
mula. 
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Considerable further theoretical work has been necessary to obtain 
the results given here. Vaulot's differential recurrence relation, which 
formulates the probability of delay at least t of a call which arrives when 
n other calls are waitmg, has no simple solution. By approximate meth- 
ods, it was possible to use a differential analyser to determine these 
jjrobabilities for small values of n. But it was not feasible in this way to 
cover the whole range of interest, and these results were supplemented 
by approximations for large n, which are described below. Finally the 
delay for an arbitrary call was obtained by summing on n* 

These results are not reported here, because the attempt to verify 
the accuracy attained led to formulation of the moments of the delay 
curves and this in turn to the representation of the curves as sums of 
exponential curves, \vith great simplification of the calculations re- 
quired. As will appear, two exponentials furnish a sufficient approxima- 
tion except for heavy traffic. 

2. DELAY CURVES 

The delay distribution on calls delayed for occupancy levels (defined 
below) from 0.1 to 0.9 in steps of 0.1 is shown in Fig. 1. The abscissae 
are derived time uiiits which seem to be natural to the problem: u = 
ci/h, with c the number of trunks and h the average holding time. The 
ordinates, on a logarithmic scale, are conditional probabilities that a 
call delayed vnW be delayed at least u, that is, values of a function 
F(u); the logarithmic scale is chosen to emphasize the dominantly ex- 
ponential character of the curves. The occupancy level a is the ratio 
a/c where a is the average call input in average holding time h. 

Fig. 1 is a master curve for all eventualities and may be changed to 
working curves for various sizes of trunk groups. For the constmction 
of these curves Table I, from which Fig. 1 was made, and which also 
compares present results ^nth those for calls served in order of arrival, 
is convenient. A more elaborate table \vill be given later. For the con- 
venience of the reader, it may be noticed that for order of arrival serv- 
ice F(u) = r""-''\ 

The striking feature of Table I is the increase in delay time for ran- 
dom service, which becomes more pronounced with decreasing F(u) 
and increasing occupancy (or traffic) level, a. The increase throughout 
the talile is an effect of the limitation to small values of F(^l). For given 

* Tlianks are due to George W. Abrams for directing this work, to Dr. Richard 
W. Hamming for transforming the equations into forms suitable for the differen- 
tial analyser and for supervising its operation, and to Miss Catherine Lennon for 
a great deal of calculation. 
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Fig. 1— Delay curves for random service. F(u) = conditional probability of 
delay at least u; u = cl/k, c = no. trunks, h = av. holding time, a = call input 
in time h, a = a/c. 

a, the delay curves for order of arrival and random service include the 
same area, which is in fact equal to the mean delay (of calls delayed), 
(^l _ a)^\ Since F(0) = 1 for both, and the random service curve 
decreases more slowly for large u, the curves must intersect at some 
point, say for it = Wo ; for w < i^ , the o.a. curve must be above the ran- 

Table I — Delay-Time and Random Service 
Delay Times, u, for given F{u) and a and for order of arrival (o.a.) and ran- 
dom service. 
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dom curve. This is shown in Fig. 2 for a. = 0.9, but the logarithmic 
scale for F{u) obscures the equality of area. 

The character of the comparison may be clearer if the picture is 
changed. Consider a department store counter with c clerks (correspond- 
ing to c trunks) in attendance. The time for a sale corresponds to the 
trunk holding time, and the rate of arrival of customers is hke that of 
call input. For service in order of arrival customers are given serially 
numbered tickets on arrival; for random service, these tickets may be 
supposed drawn from a hat, or numbered from a series of random num- 
bers, or since aggressiveness and the clerks' attention are subject to 
devious rule, it may be that no attention at all to order of service is 
equivalent to random ser^'ice. 

The fact that the average delay is independent of the order of service 
may be explained roughly by saying that the average rate at which wait- 
ing lines are removed depends only on the average rate of arrival of 
customers and the rate at which they are served. Notice however that 
service at random causes more variable delays (the second and all 
higher moments are larger than for order of arrival service). Thus with 
random service the proportion of waiting customers receiving quick 
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Fig. 2 — Comparison of delay curves for order of arrival and random service; 
a - 0.9. 
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service is increased (over order of arrival) but this is achieved at the 
cost of making other customers wait much longer. 

Service in order of arrival haS the advantage to the customer that his 
delay is independent of all who come after him, and this is particularly 
appreciated in times of heavy crowding when long delays are possible 
for random service. In Table I, these crowded conditions correspond to 
small values of Fiu) or large values of a, or both. In this picture it 
seems mtuitively clear that much longer delays are possible for random 
service, for those unlucky customers who keep missing their turn. 
(Of course, a more realistic model would also include the effects of cus- 
tomers leaving before service, a factor of considerable telephone interest 

also.) 

As noted at the start of this section, F(u) is a conditional probability, 
the probability of delay at least m of a call that is surely delayed. To 
obtam unconditional probabilities of delay, Fiu) is multiplied by the 
probability that all trunks are busy, which is the probability that a call 
is delayed. This probability is given by a well-known formula due to 
Erlang and customarily written as 



C{c, a) = 



2 c-1 

a . a . , a 



l + ^-h^+--- + 



(c - 1) ! (c - a) L 1! 2!^ ^ (c - 1)! 

+ 



(c - 1) ! (c - 



a)\ 



Tables of this function are available*. 

Finally it may be noticed here that for random service and light traf- 
fic (roughly, a less than 0.7), with sufficient approximation 

with yi = I - V^, ^2 = 1 + V^. 

* But there seems to be no extensive tabulation. However, the table for the 
Erlang B function made by Conny Palm (Stockholm, 1947) may be used with the 

relations 

1 1 1 



Cic, a) Bic, a) B{c - 1, a) 

_ a 1 — (g/c) 
c B(c, a) 

Notice that C{c, a) also has the recurrence relation 

1 -1 (c - a)(c - 1) 



C{c, a) c - 1 - a a{c - I - a)Cic - 1, a) 
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3. BASIC FORMULATION 

As noted abo\'e, the following notation is used; c is the number of 
trunks, h is the average holding time (the distribution of holding times 
is exponential) and a is the average number of calls arriving in time in- 
terval h. Then, if Fn{t) is the probability of delay at least t of a call 
arriving when n other calls are waiting, the differential recurrence 
relation given by Vaulot is 

dlM^nc _ c + a ^ a 

dt 71 + 1 ft h h 

This may be derived as follows. Consider the intei-val dt after the epoch 
of arrival of the call in question. In this interval three events may occur: 
(i) a call may arrive, (ii) a trunk may be released, or (iii) neither of 
these. The probability of a call arrival is (a/h)dt and if a call arrives 
the delay function is F„+i(i — dt). The probability of a trunk release, 
because of the assumption of exponential holding time, is (c/h)dt, and 
if a tnink is released the number of w-aiting calls is reduced by one; the 
probability that the call seizing the waiting trunk will not be the call 
in question is n/(n + 1)- Finally the probability of the third event is 
1 — (c + a)dt/h. All this is summarized in the differential relation 

F.{t) = ldt F.+,{t - dt) + -^ I dt F„.i(t - dt) 
ft n -\- 1 h 

+ ("l - ^^ dt^FM - dt) 

Passing to the hmit gives equation (1). 

Using new variables : u = ct/h, a = a/c, equation (1) may be written 
more simply as 

d_F^ ^ _n ^^^^^^^^ _ ^^ _^ ^^^^^^^^ ^ «F„+i(«) (la) 

du 71+1 

This equation is a mixed differential-difference equation of the first 
order as a differential equation and of the second order as a difference 
equation; hence three boundary relations are required. For the differen- 
tial part, it is clear that F„(0), which is the probability of some delay 
of the test call, is unity for all n in question, that is, for all integral non- 
negative n. Also F„(w) = for all negative n, is an obvious necessity, 
and, since F„ is a distribution function F„(os) = l. Finally the third 
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condition may be stated as 

lim Fn{u) = 1, all w 

The probability of delay at least u of an arbitrary call is the sum on 
n of the product of the probability that n calls are waiting when the 
call arrives and the probability, F^iu), that for this condition the call 
is delayed at least u. The first probability (for statistical equilibriuna) 
is known to be 

CI - a) C{c, a)«" 

where C{c, a), as stated above, is the probability that all trunks are 
busy; (1 — a)C(c, a) is the probability that all trunks are busy and no 
calls are waiting. Hence the probability in question, say /(w), is given by 



or by 
if 



/(„) = (1 - a)C(c,a) ^ct"FM 



f(u) = C{c, a)F{u) 
Fiu) = (1 - «) E a"^„(w) (2) 



F(w), like Fniu), is then a conditional probability, the probabiHty at 
least w of a delayed call. Notice that, consistent with this, F(0) = 1. 

It is interesting to notice that Mellor's basic equation, which in pres- 
ent notation may be written as 

'^^^ = —GM, (3) 

du n + 1 

follows from (1) if first it is supposed that F„_i(u) = F^iu) = Fn+i(u) 
and then, for clarity, Gn replaces F„ . Hence, as indicated by the third 
boundary condition, it may be expected to be useful for large values of 
n. Its solution is 

G^{u) = e-"'^"+" (4) 

A somewhat better approximation may be determined by the Mac- 
Laurin series obtained by repeated differentiation of (la) and evaluation 
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at w = 0; this is as follows 

oi {uf a{2a - 1) {uf 



FM ^ 1 - -^^ + 



(5) 



71 + 1 ' 2 (n + 1) 3! (n + 1) 

a(2a - l)(3a - 2) (u)* _ 
'^ 4! (71+ 1) ■*■ 

But this is the same* as: 

F„(«) ;^ [1 - (1 - a)u/(n + 1)]^'"^"' (5a) 

As a approaches unity, (5a) approaches (4). Equation (5a) has been 
used, for large values of a, in the direct computations mentioned above. 
It may also be noted that for a = 0, equation (la) has the solution 
(now writing F„(u, a) for F„(u)) 

F„(u,0) = *(w,n) - — — .t>iu,n - 1) (6) 

rt ~\~ J. 

where <f>(u, n) is the Poisson sum 

Finally, for completeness, note that for small values of u, the Mac- 
Lauriii series for F(u) is 

F{u) = 1 — U log r 



1 - a 

+ |a_„)[,-i^^iog^] (7) 

- ^ (1 - «) [l + 3a - (1 - a) log-i f:^] 

6 |_ 1 — a 1 n^j 

4. MOMENTS 

The A;'th moment (about the origin) of the delay density function 
which is —F'{ti) (F(u) itself is a distribution function is defined as 



M* = / u'[~F'(u)]du, 
Jo 

= k I u''~^F(u) du, k > 0, 

Jq 



(8) 



the last by integration by parts. 

* G. W. Abrams is due credit for noticing this. 



(10) 
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Following (2), this may also be written as 

so 

M, = {I - a) ^ ci"mn,k , (9) 



with 

m„.ft = [ u[-F'M] du, 

= k I vl'-'FM du, k> 0. 
Jo 

First, notice that 

m„.o = -f F'Mdu = FM = I.- 
Jo 

hence 

so 

j¥o = (1 - «) E «" = 1. 

showing that F{u) is properly normahzed. 

Next, by integrating both sides of (la) with respect to u from to m , 
and using the second form of (10) (with A; = 1) 

-(n + 1) = ?tm„_i.i - {n + 1)(1 + a)m„,i + (n + l)am„+i,x (H) 

In the same way, after first multiplying (la) throughout by u ~ , it is 
found that 

-fc(n+ l)m„,t_i ■ 

= nm„_i,i - (n + 1)(1 + a)mn.k + (" + l)am„+],A 

Unfortunately, neither (11) nor any other instances of (12) have simple 
solutions; nevertheless they may be used to determine Mk. 

Consider first the simplest case. Mi . If (11) is multiplied throughout 
by a" and summed on n, the result may be written 

— i/io = aLn — (1 + o:)-^ii + -^n ~ -^oi 

= *~-^i 

where for convenience in writing and of later notation 

Loi = E ""mn.i = (1 - a)-' Ml 
Lii = Y^{n -\- l)Q!"m„,i 

u = E (« +!)«" = -D E """"^ = (1 - «)"' 



(13) 
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and D = d/da. Hence 

Ml = (1 - a)"' 

This is the mean delay of calls delayed and as mentioned above is the 
same as for service in order of arrival. 

In the general case*, the following notation is convenient 

Lok = J2 <x''mn,k = (1 - ay^Mk 
Ljk = 'Z(n-\- l)(n + 2) ■ ■ ■ (n + j>"m„,fc 
Using the relations 

n(n + 2) • ■ • (n + j) 

=-n(n+l) ■•• (n+i-l)+(y-l)n(n+l) • • . (n -\- j - 2) 
+ • • • + (i - l),n(n + 1) . . . (n + _^- - i - 1) -f . . . 

+ 0' - iy.n, 

n(n -\- 1) ■•■ (n + j - I) 

= (n + l)(n + 2) • • • (« + j) - j(n + 1)(^ + 2)(n + j- 1) 
mth 

0"- 1).'= (J - 1)0' -2) ■••(j-i), 
the summing of (12) is found to result in 

kLj,k-i = [j - (j - l)a]Lj-i,k - oi[(j - lhLj-2.k , ^ 

(14) 
+ (j - l).Lj-,,k + ■■■ + (>- l).Lj-i.k +■■■ + (;- ly.Li.k] 

But this may be simplified by multiplying tlirough by j and sub- 
tracting from the same equation \\ith j replaced by j + 1; the result is 

(j + 1 - jo^Ljk - fLj-i,k = kLj+i,k~\ — jkLj,k-i (15) 

Notice that for j = 0, k = 1, Lm = Lio , as in (13). Notice also that 

<^^%-o = E (« + 1) ■ • ■ (n + j>"+''-^ 

= Dj^{n-\-l)--- (n-\-j ~ l)a"+^' 

= DCa'Ly_i,o) 

so that 

Ljo = iL,_i,o + «Z)L,-_,,o = ;!(1 - ay'-' 
* This procedure is the development of a suggestion made by S. O. Rice. 
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Then the ratio 

Lo,{L,,r' = (1 - ^xf'-'U/kl 

= (1 - afMJk\ = R, 

is the ratio of these moments to those for order of arrival service; the 
last relation is a definition. In the same way the ratio 

might be considered, but to avoid fractions the following somewhat odd 
change of variables seems convenient: 

{'V')9kLik = psf^Lj^f..^ _ (17) 

where ga = gi = I and 

ff2fc = (2 - af-\S - 2af'-' ■■• (k+l-Jca) 
g,,^, = (2 - af(3 - 2af-' • ■ ■ (fc + 1 ~ kccf 

Notice that 

gu+iig^k)-' = ff2.(s2.-i)-' = (2 - a)(3 - 2a) ■ ■ • (fc + 1 - M = D^, 

the last being a definition, again. 
Since 

(1 — a)Lic,o = kLk-i.o 
it follows from (15) that 

(i + 1 - J«)pyt - y(i - «)pj-i.t ,jg^ 

= (gk/9t-i)[U + i)p>+i.*-i - J'd - «)p/,fr-i] 

By taking differences of this equation and writing 
go* = Pot 

Qik = pit — Pok = Apofc 

q2k = p2t — 2pift + pofc = A pot *^ 

gjk = Aqj-i,k = A^pofc 
a somewhat simpler recurrence relation is found to be as follows 

(i + 1 - i«)g;fc = igk/gk^i) , 

[jaqj-i,k-i + (i + 1 + JodQj.k-i + 0' + l)5>+i.t-i] 
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Since pjo = 1, ally, gw = 1, and gyo = 0, j 9^ 0. From these boundary 
conditions, it follows at once from (19) that 

gsk = 0, j > k. 

By comparison of (17) and (16) 

gkRk = pok = qok 
A short table of the g's is as follows: 



J/* 





I 


2 


3 




1 

2 
3 


1 





1 

a(2 - a)-i 




2 

4a (2 - a)-i 
2a»(3 - 2a)-i 



2{2 + a) 
2a (18 - 5« 
2«ni8 - 7« 


- 2a^)(3 - 


- 2a)-* 



Continuation of this leads to the values of Rk listed in Table II. Notice 
that for «= 1, by (18) 

P^k (1) = U + l)p/+i,.-i (1) 

= (j + i)(i + 2)pj^,,k-2 (1) 

= U + 1) (j + 2) ■ ■ • (i + k) 
since fft = 1 for a = 1 and pyo = 1 , all j. From this 
PokH) = gMRk{l) = Rk(\) = k\ 
On the other hand, for a = 0, gjk = 0, j > and, by (19) 

gok(0)/gk(0) = ?o.t-i(0)/ff,_i(0) 
so that 

Rk{0) = i2*-i(0) = Rt{0) = 1 

5. MELLOR APPROXIMATION 

It is useful to have the moments of the distribution corresponding to 
the Mellor approximation, since they serve as a guide. Here, following 
equation (4) 



Fiu) = (1 - a) i: a 



n -ii(n+l)-l 



(20) 
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and 

exp - xM = E Mk{~xY/k\ 



= (!-«) r du e^^" i; (71 + l)-^a>.-"<"+"~^ (21) 
Jo 

= (1 - «) f; «"[i + x{n + i)r 

D 

Hence 

jlf, = fc!(l - a) i; (n+ 1)V (22) 



These moments are expressible in terms of polynomials associated with 
the distribution of permutations into classes according to the number of 
readings left to right necessary to find the elements in standard order. 
Indeed the ratio 

r,(a) = M,a- af/h\ 

has the recurrence relation 

r,+i(«) = {ka + l)r,(a) + a(l - a)r'fc(a) (23) 

and the first few values are as follows 

ri = 1 rs = 1 + 4a 4- a 

r^ = 1 + « n = 1 + lla + 11a' + a 

r^ = 1 + 26a + 66a' + 26a^ + a 

Notice that n{Q) = 1, ^{1) = /:!, just as for the precise results. 

6. EXPONENTIAL SUMS 

The shape of the delay curves, from direct calculation, and also from 
Mellor's results, suggests representation in exponential sums. If 

F{u) = Ai e-"-"""^' + A2 e-'^-"'"'-+ ... (24) 

then 

M, y^^^' = A,x\ + A^x^ + ■ - . (25) 

by a simple calculation. For h exponentials, 2fc moments (including 
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Mo) may be fitted exactly by solution of 2k equations of form (25), as 
will be shown. 

The first approximation {k = 1) is the order of arrival curve, say 

which has Ai = Xi — 1, Ak = Zk = 0, k > 1, and matches Mo and Mi . 
The next approximation (fc = 2) is determined from equations 

Ai + .4, =1 

AiXi + A2X2 = 1 

Aix\ + A2XI = Rz 

Aix\ + A^xl = R3 

Eliminating A2 from successive pairs, 

Ai (xi - 0:2) = 1 - X2 

AiXi (xi — .1:2) = R2 — X2 

Aixl(xi — xt) = Ri — R2X2 

Eliminating Ai from these, 

Xi + X2 — X1X2 = R2 

(xi + X2)R2 — X1X2 = Ri 



(26) 



or, writing ai = Xi -{- X2 , 02 = XiXa , so that x — aix + 03 = (x — Xi) 

(x - X2) 

fll — (32 == R2 

(26a) 

CI1R2 — (l2 = /l3 

From the first of the second set of equations, and from symmetry (or 
from Ai -\- A2 = 1) 



A 1 — X2 

Ax = 

Xi — X2 
A2 = 

X-> — Xi 

Taking R2 and Ri from Table II, it turns out that 
Xi~' = 1 - V^ = 2Ai 
X2~' = 1 + VW^ = 2^2 



(27) 



(28) 



.t*^ 
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Table II - Moment Ratios, Calls Served at Random 



Ri = 1 

2 



R,= 
R, = 
R,= 
R,= 



2 - a 

2(2 + a) 

(2 - «)" 

4(6 + 5a - 4a' - a^) 

(2 - a)K3 - 2a) 

4(36 + GQq - 59g' - 24a^ + 15a* + 2«') 

(2 - «)*(3 - 2«)» 

8/a(a) 



(2 - a)''(3 - 2a) =(4 - 3«) 

8/7(«) 



(2 - aH3 - 2a)H4 - 3a)» 
/.(a) = 432 + 972fl - 2016a« - 437a» + 1790«* - 528a» - 196a« + G7a' + 6a« 
/7(«) = 10368 + 34560a - 89208a" - 32772a' + 177926a^ - 104287a* - 29260a<' 

+ 43876aT - 9158a» - 20390* + 588«i« + 36ai 



and the second approximation is 

2¥^{u) = (1 - Va/2) e ^ 

+ (1 + VW2) .-""-"'"+^'^'=' 

which turns out to be a good fit for a roughly less than 0.7. Curiously 
the corresponding Mellor approximation has a more complicated ex- 
pression. 

Following the same procedure for three exponentials, it turns out 
that the correspondent to the set of equations (26a) is 

aiK3 - 12^2 + as = R^ (30) 

with ax = Xi +. Xi-V xz,a^ = 0:1X2 + x^Xz + x^x^ , 03 = xia:2X3 , that is, 
the symmetric functions. 
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Using Table II for values of the R's, it is found that 
ai = (18 - 7a - 2a)(2 - a)'\3 - 2a)~^ 
02 =18 (2 - a)"\3 - 2a)-' (31) 

ag = 6 (2 - a)~\S - 2ay^ 

X\ , xz and Xz are then the roots of the cubic equation 
X — aix -\- ttiX — as = 
The coefficients Ai , i = 1, 2, 3 are determined from equations like 

A> = ^^ - (^^ +/^) + ^;^^ (32) 

[xi — 0:2) (a;i — X3) 

For the fourth approximation, matching 8 moments, the equations 
for the symmetric functions are 

O1B3 — diRi ~l~ fl3 — fli = Ri 

UiRi — aafla + diRi — O-i = R& 

(33) 
ttiRi — ttzRi + a-jRi — a\Ri = Ra 

ttiR^ — OiRf, -\- ttsRi — a^Rz = Ri 

and Xi , X2 , 3:3 and x^ are roots of the quartic equation 

X* — aix^ + a2X^ — ttao; + 04 = 

Coefficients Ai are determined from equations like 

, ^ Rz — (xb + a:a + 3:4)^2 + (3:23:3 + a^aa:^ + 3^3X4) ~ 3:2X3X4 ,^^-. 

(Xi — X2)(xi — X3)(xi — Xi) 

It may be noted that 

3:2 + X3 + X4 -= Oi - Xi 

X2XS + 3:2X4 + X3X4 = 02 — Xi(ai — Xi) 

X2X3X4 = 03 — Xi[a2 -- Xi(a — Xi)] = 04X1' 

which gives the general structure. 

It is worth noting that equations (33) may be used to determine the 
R's if the a's may be determined otherwise. As a matter of fact, they 
have led to the determination of R^ and /^y in the following way. The 
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results for k = 2 and 3 suggest that 

a* = 4!(2 - a)"'(3 - 2ar'(4 - 3a)"' 
az = 'itti 
Then by the first two of equations (33) 

aiBz — aiRi = Ri — a-i -{- at 
aiRi — a-^i = Rh — a-iRi + a* 
the solutions of which are 

ai = 4(24 - 23« + Za') [{2 ~ a)(3 - 2«)(4 - Sa)! 

02 



2(72 -- 23a - 10a' - 3a')[(2 - a)(3 - 2a)(4 - 3a)] ' 



By the last two of equations (33), R^ and Ri are determined to be the 
vahies given in Table XI, which have been verified independently. Note 
that for a = 0, both R^ and Ri are 1, and for a = 1, 2^6 = 6!, Ri = 7! 
Table III tabulates, for k = 2 to 5, for convenience in avoiding frac- 
tions the symmetric functions hk, related to those above by 

hki = Dhttj 

with, as before, 

Dk = (2 - a)(3 - 2a) ■ - - [k - (k - l)a] 
and oo = 1- The functions for A; = 5 were obtained by a process like 

Table III — Symmetric Functions for Exponentl\l 
Sums of Calls Served at Random 



it = 2 &20 = 2 - « -t - 3 6„ = 6 - 7a + 2a« 

6„ =4 fta, = 18 - 7« - 2«' 

b22 = 2 bii = 18 

6aj = 6 

it = 4 bio = 24 - 46a + 29«2 - 6a' 

6„ = 96 -- 92q + 12a= 

bt- = 144 - 46q - 20a= - 6a= 

6(3 = 96 

644 = 24 

k = 5 bfio = 120 - 326a + 329a2 -- 146a' + 24a< 
b» = 60U - 978a + 329a^ + 146a» - 72a* 
bi2 = 1200 - 978a - 172a^ + 78a» + 72a^ 
661 = 1200 - 326a - 172a2 - 78a' - 24a» 
bst = 600 
6sfi = 120 
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Table IV — Symmetric Functions for Exponential Sums, Mellor 
' Approximation 

k = 2 fli = 3+a fc = 3 oi = 6 + 3a 

ai = 2 Oj =: 11 + 5(^ _|_ 2aa 

aj = 6 

A- = 4 fli = 10 + 6a 

aa = 35 + 26a + lla^ 

fla = 50 + 26a + 14a^ + 6«' 

a, = 24 

/. - 5 ai = 15 + 10a 

02= 85 + 80a + 35a2 

fla = 225 + 200a + 125a^ + 50a^ 

at = 274 + 154a + 94a^ + 54a» + 24a* 

ao = 120 

fc = 6 oi = 21 + 15a 

fla = 175 + 190« + 85a= 

as = 735 + 855a + 585a2 + 225a3 

04 = 1624 + 1604a + 1194a" + 704a« + 274a' 

at = 1764 + 1044a + 684a' + 444a' + 264a* + 120aS 

a& = 720 



thiit sketched above, and without determining Rg and Rg . Notice that 
hj = k\{)), a = 

= JK))\ « = 1 . 

which may be proved independently. All values in Table III satisfy 
the retiurrence relation 

hj = [k - {k - l)a]b,-,,,- + [k -\- {k - l)a]6*_i,y_i (35) 

— (k ~ 1) abk-2.j~2 

which also satisfies the boundary relations for a = and 1 given above 
for all values of k. 

The corresponding symmetric functions for the Mellor approximation 
are given in Table IV. These have the recurrence relation 

flfcj = ak-ij -\-[k-\- (k — l)a]ak-i.3-i — {k ~ l)W-2.j-2 (36) 

For a ^ 0, the values are the signless Stirling numbers of the first 
kind, that is, the numbers given by the expansion of 

(1 +a-)(l + 2x) ■-■ (1 -\-kx). 

For a = 1, the results are the same as for the exact ease, as given above. 
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Table V — Approximations to Delay Function F{u) for Random 

Service 



Two Exponentials 



0.1 
0.2 
0.3 
0.4 
0.6 
0.6 
0.7 
0.8 
0.9 



.3590 
.3490 
.3392 
.3292 
.3190 
.3085 
,2978 
.2868 
.2756 



.1351 
.1344 
.1332 
.1315 
.1293 
.1265 
.1232 
.1193 
.1148 



.0220 
.0256 
.0291 
.0325 
.0357 
.0386 
.0412 
.0434 
.0451 



.0041 
.0058 
.0079 
.0101 
.0125 
.0151 
.0177 
.0203 
.0229 



.0008 
.0015 
.0023 
.0033 
.0046 
.0061 
.0078 
.0097 
.0118 



.0002 
.0004 
.0007 
.0011 
.0017 
.0025 
.0035 
.0047 
.0061 



.0001 
.0002 
.0004 
.0006 
.0010 
.0015 
.0022 
.0031 



.0001 
.0002 
.0003 
.0004 
.0007 
.0011 
.0016 



Three Exponentials 



0.1 


.3586 


.1354 


.0219 


.0040 


.0008 


.0002 






0.2 


.3491 


.1356 


.0254 


.0057 


.0014 


.0004 


.0001 




0.3 


.3393 


.1358 


.0288 


.0074 


.0022 


.0007 


.0002 


.0001 


0.4 


.3291 


.1360 


.0322 


.0092 


.0030 


.0011 


.0004 


.0002 


0.5 


.3186 


. 1363 


.0358 


.0112 


.0040 


.0016 


.0007 


.0003 


0.6 


.3071 


,1359 


.0392 


.0133 


.0050 


.0022 


.0010 


.0005 


0.7 


.2951 


.1354 


.0428 


.0156 


.0063 


.0028 


.0014 


.0007 


0.8 


.2822 


.1344 


.0466 


.0181 


.0077 


.0036 


.0018 


.0010 


0.9 


.2683 


.1325 


.0504 


.0210 


.0094 


.0045 


.0023 


.0013 



7. numerical results 

Table V gives both two-exponential and three-exponential 4 decimal 
approximations to the delay function F{u) for 



and for 



a = 0.1(0.1)0.9(0.1 to 0.9 in steps of 0.1) 



«(1 -a) = 1(1)2(2)14, 



in the same abbreviated notation.* The variable y = u(l — a) is in- 
troduced to reduce the spread of these tables. It will be noticed that, 
as expected, the two orders of approxunation agree closely for small 
values of a; indeed, only for the three largest values of a are the dif- 
ferences appreciable from the engineering standpoint. 

* The results for two exponentials, some of those for three-exponentials, and 
all special results given below, have been obtained bv Miss Marian DarviUe, 
whom I also thank for her careful drawing of the curves. The entire three -exponen- 
tial table has been computed independently by Miss Lennon. 
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For a = 0.9, results for four exponentials have also been obtained 
and compare with those of Table V as follows (k = number of ex- 
ponentials) : 



k 


V 


1 


2 


4 


6 


s 


10 


12 


14 


2 
3 

4 


.2756 
.2683 
.2748 


.1148 
.1325 
.1402 


.0451 
.0504 
.0483 


.0229 
.0210 
.0195 


.0118 
.0094 
.0091 


.0061 
.0045 
.0047 


.0037 
.0023 
.0026 


.0016 
.0013 
.0015 



It is somewhat surprising that two exponentials should do as well as 
they do for large values of v (in fact for v = 12 and 14 better than three) ; 
a similar behavior appears in the following comparison of approxima- 
tions on the Mellor basis, again for a = 0.9 



k 


V 


1 


2 


4 


6 


8 


10 


12 


14 


2 
4 

6 


.2725 
.2671 

.2777 


.1115 
.1379 
.1408 


.0446 
.0502 

.0477 


.0237 
.0207 
.0205 


.0129 
.0097 
.0102 


.0070 
.0051 
.0054 


.0038 
.0029 
.0031 


.0021 
.0018 

.0018 



From these comparisons, it appears a relatively small number of 
exponentials is sufficient for engineering purposes. The curves of Fig. 1 
are those for three exponentials, for uniformity. 
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